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Abstract 

In this paper, the effects of boundary and connectivity on ideal gases in 
two-dimensional confined space and three-dimensional tubes are discussed in 
detail based on the analytical result. The implication of such effects on the 
mesoscopic system is also revealed. 
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1. Introduction 

With the increasing interest in studying the quantum effects in mesoscopic systems 

[1, 2] which are so small that the boundary effect gets important and can not be 

neglected any more, the question of how to solve the sum of the states precisely, 

or, how to take the boundary effect of such system into consideration, is raised 

naturally. The properties of some systems found recently are shape dependent and 

sensitive to the topology [3, 4, 5]. 
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The key problem in statistical mechanics is to solve the sum over all possible 
states. In principle, when an ideal gas is confined in finite volume, the spectrum 
of single-particle states will be determined by the configuration of the boundary. 
In case the boundary is irregular, it is impossible to get the sum over all possible 
states exactly. However, if the mean thermal wavelength of the particles is much 
shorter than the size of the system, as an approximation, one can assume that the 
spectrum of single-particle states be continuous while the total number of states be 
independent of the shape of the boundary and simply proportional to the volume 
(in three dimensional case) or the area (in two dimensional case) of the system. 
In other words, if the thermal wavelength of particles is very short in relative to 
the size of the volume or area which the system occupies, the effect of boundary 
configuration on the spectrum can be ignored. Historically, such an assumption was 
advocated by radiation theory of Rayleigh- Jeans. It also aroused great interest to 
the mathematicians, and finally was proved by Weyl mathematically [6]. 

The above problem in statistical mechanics is related to such an inverse problem 
in mathematics, that is, if it is possible to determine the metric and topological 
features by the knowledge of the spectrum, which is still an open question. 

After Weyl's leading work[6], some progress has been made in mathematics [7]. 
Now we know that, in two dimensions, for the eigenvalue problem 

]-V 2 U + nU = in ft with U = on T, (1) 

where Q is the region bounded by curve T, the total number of eigenstates N can 
be written as [7] 

Here f2 is the area of Ct, L the length of T, and r the number of holes in f2, while 
{Hn\ is the spectrum of eigenvalues of the system. In the short-wavelength limit, 
Weyl proved for two-dimensional case that iV ~ i.e., the total number of states 
is proportional to the area of the region, which is a fundamental assumption in 
statistical physics. 
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Nevertheless, Weyl's result omits the last two terms in eq.(2) coming from the 
perimeter and the number of holes and hence loses part of the information related 
to the geometry of the region. Here, the relation between the perimeter L and the 
area Q contains some information of the shape of the region while the number of 
holes r reflects the connectivity. 

In this paper, the effects of the boundary and the connectivity on the quantum 
statistics are discussed. In Sec. II, the analytical result with the geometry effects 
being considered for ideal Bose and Fermi gases in two-dimensional space is given, 
based on which the relevant physical indications of such effects are revealed. In 
Sec. Ill, the quantum statistics in there-dimensional tubes is discussed. The con- 
clusions are summarized in Sec. IV while some expressions of useful thermodynamic 
quantities are given in Appendix. 

2. Statistics in two-dimensional space 

Obviously, eq.(l) is just the Schrodinger equation for free particles in a two-dimensional 
container. Based on eq.(2), we have 

where A = h/^/2nmkT = nm is the mean thermal wavelength, e s the 

energy of a free particle, and the subscript s labels different states. In macroscopic 
systems, the area (volume) of the system is usually large enough so that the influence 
of boundary can be ignored. However, it may not be true for some special cases. 
The above formula is a more precise approximation which includes the influence of 
the shape and the connectivity. It means that the geometry of the container may 
result in observable effects in physics. The effects described by the last two terms 
of eq.(3) can be expected to be observable in such cases that, 1) the area (volume) 
of the system is small, 2) the area (volume) is not small but the boundary of the 
system is so complicated that L ^> \/VL or the system is in multiply connected space. 
For ideal Bose and Fermi gases in two-dimensional space, when the influence of 
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boundary and connectivity is considered, the grand potential of the system reads as 
[8] 



\nE = Tj2 ln ^Tze~^ s ). (4) 

s 

In this equation and following, the upper sign stands for bosons and the lower sign 
for fermions. Substitute eq.(3) into eq.(4) and expand it as a series of ze^ 1363 , we 
have 

in2 = EE[(±ir +1 ^(^- fe ) n ]- (5) 

s n U 

Of course, such expansion is valid only for < z < 1. 

By using the eq.(3), we can perform the summation over s: 

Now, the sum over n gives Bose-Einstein or Fermi-Dirac integral 

Here, we introduce a function h a (z) which equals to Bose-Einstein integral g a {z) or 
Fermi-Dirac integral f a (z) in Bose or Fermi case, respectively. Finally, the grand 
potential can be expressed as 

InS = ^h 2 (z) - \jhs_(z) + (8) 

In comparison with the grand potential lnS/ ree = Qh 2 (z) for ideal gas in two- 
dimensional free space without boundary, grand potential eq.(8) depends not only 
on the area but also on the boundary as well as the connectivity of the space. The 
term which is proportional to the perimeter L of the area, represents the influence 
of the boundary while the another term reflects the effect of the connectivity. 

In addition, it is easy to see from eq.(8) that, grand potential of a system in 
confined space is less than that in free space since the signs of these two terms are 
negative (when r > 1). It means that the existence of a boundary and holes tends 



4 



to reduce the number of states of the system. This is just because in free space the 
energy spectrum is continuous while in confined space the spectrum gets discrete. 
Thus the number of modes in confined space is less than that in free space. 

The expansion in eq.(5) requires < z < 1. In Bose-Einstein statistics, such a 
constraint on the fugacity z is naturally satisfied. In Fermi-Dirac statistics, we have 
< z < oo. Strictly speaking, the grand potential InS can not be expanded as a 
series of ze~^ s when z > 1. However, the first term of InS in eq.(8) is just the 
grand potential in free space lnSj ree though the expansion is not rigorous. Since 
such a treatment provides a correct result in free space, we may expect that the rest 
two terms of InS, which describe the contributions of boundary and connectivity, 
are also valid for z > 1. 

When the thermal wavelength A is much shorter than the size of container, 
particle can not feel the shape of boundary. In contrast, for low frequency waves 
whose wavelengths are comparable to the size of container, its spectrum will seriously 
depend on the boundary. Therefore, when the ratio between the wavelength and 
the size of container is extremely small, a good approximation can be given by the 
first term of eq.(8) which is just the result proven by Weyl, and is only related to 
the area of the container. However, when the ratio is not negligible, the influence 
of the container geometry has to be considered, which is provided by the last two 
terms in eq.(8). 

One point must be emphasized here. When solving statistics problem, one has 
got used to assume that, the number of states be only proportional to the area 
and neglect the effect of boundary. Eq.(8) provides a more complete and precise 
approximation. 

Following general procedures, the relevant thermodynamic quantities of two- 
dimensional ideal gas can be achieved easily although the extra two terms in eq. (8) 
make the derivation get tedious. Eliminating z in the two equations 



we obtain the equation of state of the ideal gas. The other thermodynamic quantities 
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1 — r 
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are given in Appendix. 

The effect of boundary on the thermodynamic quantities (see Appendix) is of the 
order of L/(y/Ny/U). The factor L/VU reflects some information of the shape of 
the two-dimensional container. If the container shape is close to circle or square, the 
ratio L/y/Ti is of order 1. Otherwise, for example, when the shape of the container 
is very complex, L/y/Q can be large and hence the boundary effect will become 
significant. 

Besides, the boundary effect is suppressed by the factor l/y/N. In macroscopic 
systems, the contribution of boundary is strongly suppressed. But, the suppression 
may not be so serious for mesoscopic systems, and so that is expectable to observe 
the boundary effect in such systems. 

The influence of connectivity is of order (1 — r)/N, so one may observe the effect 
of connectivity when r is comparable to N. This result implies that in some porous 
media such effects may not be ignored. 



3. Statistics in there-dimensional tubes 

Base on eq.(3), we can also calculate the boundary effect of a three-dimensional 
ideal gas in a long tube, of which all transverse cross sections keep the same. 

The ^-component of the momentum p z is continuous since the length of the tube 
L z is made sufficiently large. This allows us to convert the summation over p z into 
an integral. Then, we only need to perform the summations over p x and p y , which 
can be achieved by using the same procedure above. The grand potential is 



InS = £(±1)" +1 — £ e - n/3£s 

n n s 

= E(±l) n+1 - /^V^i £ e~^. (10) 
Working out the summations over p x and p y by use of eq.(3), we obtain 

lnE=M ¥2 )-iM fc2(2) + i^ ¥2 ). (u) 
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where Q and L denote the area and perimeter of the transverse cross section of the 

tube, respectively. 

Directly, we can obtain the equation of state 

PV L Z Q 1L Z L 1-rL 

{ N = ^-h,(z)--- jr h 1 (z) + ——h h (z). 
The thermodynamics quantities are listed in Appendix. The Bose-Einstein conden- 
sation in a three-dimensional tube will be discussed in detail elsewhere. 

A similar analysis indicates that the effects of boundary and connectivity in 
three-dimensional tube are of the same order as those in two-dimensional space. 



4. Conclusions and discussions 

In conclusion, the effects of boundary and connectivity on the statistical mechanics 
of ideal gases in two-dimensional confined space and in three-dimensional tubes are 
discussed. 

In ideal gas theory, one replace the summation over states by an integral over 
momentum: J2 S — > ^/^ri where d denotes the dimension. This replacement is 
based on the assumption that the momentum is independent of the boundary of 
container, which is valid only when V — > oo, i.e., there is no boundary. In other 
words, such a replacement is equivalent to the assumption that the number of states 
is proportional to the area (in two dimensions) or volume (in three dimensions) of 
the container. In fact, this treatment is an approximation that ignores the influence 
of boundary. 

Our analysis shows that, the influence of boundary and connectivity is of the 
order L/(y/Ny/Ti) and (1 — r)/N, respectively. In many cases, such influences are 
of order 1/y/N and l/N, therefore, negligible. However, when these factors L/\/~VL 
and r are much bigger than 1, the effects may be observable. The factor L/VU 3> 1 
corresponds to such cases, for example, the boundary of the region is very complex 
or the two-dimensional container is long and narrow, r ^> 1 corresponds to the 
case that there are many holes in the region. In the meanwhile, the suppression by 
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l/y/N or 1/N could also be reduced in mesoscopic systems, so we can expect that 
the effect of container geometry may be observable in mesoscopic scale, especially 
in containers with complex boundary or in porous media. 
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Appendix: Thermodynamic quantities 

1. Thermodynamic quantities in two dimensions 

Internal energy: 

U h 2 (z) 1 L h 3/2 (z) i /2 

ff 2 ~ ~ 7= T^TTTT, 7 a 2 ■ 



NkT hiz) 8VNVnh{ /2 (z) 
Free energy: 

F , fh 2 (z) 1 L h 3/2 (z) 1/2 l-r : 
- \nz - — -d 2 - —7=^= ' <r 2 + -zrrh^z) 



NkT yh^z) A^Vnh\ /2 (z) QN 

Entropy: 



S n h 2 (z) 3 L h 3/2 (z) 1/2 1 -r 

-777" = 2- — —cr 2 — mz = — (To H 77-/11(2) 

iVA; hiz) 8^V^hl /2 (z) 6N n } 



Specific heat: 



where 



Cy 

Nk 




1 h{ /2 (z)h l/2 (z) 

V2- 



cr 2 



1 ~~ liW^o(z) 



1 4- -J-H. h ll2fl _ l-r I, 

1 ^ 64JV n 6JV ""Ol^J 



1 L h i/2(z) 
SVWVU h\ ,2 (z) ) 
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The following relation is used to calculate Cy 



dz z hi(z) 
df = ~Th {z) V2 ' 

2. Thermodynamic quantities in three dimensions 

Internal energy: 

U 3 /i 5/2 (z) 1 L\I*L h 2 {z) 2/3 l-r L 2 J 3 2/3 1/3 

NkT 2h 3/2 (zf 3 47VV3 Q2/3 fgfcf 3 t 12iV 2/3nl/3 n 3/2l^3 ■ 

Free energy: 

i 71 , 1 1 L l J ? 'L h 2 (z) 2/3 

lnZ ~U /^ ff 3 - TTnTa n2 / 3 ,2/3/^3 



NkT l h 3/2 ( Z ) 6 4iW3 fi2/3 ^2/3 (z) 



1 l-r L 2 / 3 2/3 1/3 
+ A^-^^ ft 3/2l^ ff 3 J- 



Entropy: 



S _5h 5/2 (z) 1 LV3 L ^) 1-rLf 2/3 1/3 

Nk~2h 3/2 (Z) 3 27VV3 Q2/3 ^2/3^^3 +4^2/3^1/3^3/2^3 • 



Specific heat: 



r» /\r/r ^ h 5/2 {z) 9 ^3/2(2) 

C y = iv/c{a 3 [— - - -773- 



4 h 3/2 (z) 4'°h 1/2 (zy 
1 Lltj* 2/3r h 2 (z) 3 



where 



03 



7VV3 fi 2/3 ff 3 L^/^^ 8 ^3 
+ V"iV^T^3-°" 3l/3[ I (1 " ^^3/2(^)1} 



- f 1 _ ^jzlM^Mj , ^ (l-r) 3 h? /2 (z) \ ' ,-3 
V ^ 72N h 3/2 (z) n 291QN 2 h 5/2 (z) ) Kx ' 



9 



I _ 1 L^J'^L hjjz) 1 i. — i u - ■ ■ .1- / - ' - ' _i_ 

6JVV3 h 2/3 (2) ^73 + 18JV2/3nl/3 fc l/3 w T7J 



1-r L^ /3 fel/2(g) 1 



^3 



1 



l-r L^ 3 1 



J. \" 1 ! < >- :1 ^l /2 (z) a 1 / 3 + WV273 Ql73 ft 1/2 (2) CT 2 / 3 



1 1 L^L 1/3 l-r I^/ 3 hi/ 2 (z) 1 



612iVV3 fi 2/3 ^2/3 (z) 



1 1 

2 + 2 



1 



1 + 



1 L Z L 5 h\{z) 
432iV VI 2 h 2 3/2 (z 



187VV3Q1/3 
1/3 



£3 



^4 
£ 5 



We also have 



_ l-rL z L h^hy^z) (l-rf L 2 z h\ /2 {z) 
72N Q h 3/2 (z) 291QN 2 VI h 3/2 (z)' 



1 - 



L z hj /2 (z) 



27 N L fn(z) 



dz 
df 



3 z h 3/2 (z) 
'2Th 1/2 (z) 
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